function ODfilt = FilterImage(varargin)

if nargin == 1
    ODraw = varargin{1};
    N = 1;  %amplitude of 2D gaussian filter
    sig = 10;  %size of 2D gaussian filter
elseif nargin == 3
    ODraw = varargin{1};
    N = varargin{2};
    sig = varargin{3};
else
    error('Number of arguments should be 1 or 3')
end

% Filter image to find ROI
ODfilt = conv2(ODraw,gaussian2d(N,sig),'same'); % Low-pass filter image
end
